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Recent evidence for a CP violating asymmetry in the scmilcptonic decays of B s mesons cannot 
be accommodated within the Standard Model. Such an asymmetry can be explained by new physics 
contributions to AB = 2 components of cither the mass matrix or the decay matrix. We show that 
mixing with a hidden pseudoscalar meson with a mass around 5 GeV can result in a new CP 

^\ violating contribution to the mixing and can resolve several anomalies in this system including the 

width difference, the average width and the charge asymmetry. We also discuss the effects of the 
hidden meson on other b physics observables, and present viable decay modes for the hidden meson. 

f — We make predictions for new decay channels of B hadrons, which can be tested at the Tevatron, 

< ^> the LHC and B-factories. 
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1 Introduction 

In recent years, we have seen many new experimental measurements in the neutral B meson system 
both from the B-factories and the Tevatron. While most observations agree well with the Standard 
Model (SM), there are a few disagreements at the 2 or 3a level. Recently, the D0 collaboration has 
announced evidence for a charge asymmetry in the number of like-sign dimuon events pQ , which can 
be interpreted as a CP violating asymmetry in B s meson oscillation rates and semileptonic decays. 
Such an asymmetry could result from a phase difference <pf between r 12 and m 12 , 

^ = arg(-mf/lf), (1) 

where r 12 is the off diagonal term in the B s decay matrix resulting from interference between B s and 
B s decays, and m 12 is the AB = 2 mass mixing term. In the SM, the dominant contribution to ml 2 is 
proportional to the weak Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix elements (V tb V^ s ) 2 
while V 12 is approximately proportional to (V cb V c * s + V^) 2 - [8m 2 c /(3m 2 )](V cb V c * s + V^V^V^V*. 
CKM unitarity constrains V ts V t * b to equal -V cs V* b - V us V* b , predicting 0[3] 

^i(SM) « ^| x /3 S W(SM) = f^f x ar § (-TrST) = °- 0042 ± °- 0014 > ( 2 ) 

^ ra b Am b V V cs V cbJ 

therefore the SM prediction for the asymmetry is unmeasurably small, making the charge asymmetry 
an interesting place to look for new physics. 

Several other observables in the B s system are in marginal disagreement with the SM. From the 
particle data group (PDG) [3], the average lifetime of the neutral B s mesons is tb s = 1.472^ 0026 ps 
and the lifetime of the neutral B^ mesons is r^ d = 1.525 ± 0.009 ps. The ratio of those two lifetimes is 
r B Jr Bd = 0.965 ± 0.017, which exhibits a 1.8a deviation from the SM prediction of 1.00 ± 0.01 [6]. 
Using the measured width of the B& (we assume that new physics does not modify the mixings in the 
Bd system throughout this paper) , the SM model prediction of the average width of B s is 

f s (SM) = 0.654 ±0.008 ps" 1 , (3) 

while the measured average width is T s = 0.680 ±0.012 ps -1 . 

The time-dependent CP asymmetry in B s — > J /ip $ decay also determines various mixing parame- 
ters in the B s mesons. Using the combined results from D0 and CDF [3|8j with 0-2.8 fb _1 luminosity, 
the two extracted quantities are Ar s = 0.154±°;[j^ ps" 1 and /3 S W = 0.39±g;i|, where 2/3/ /,/ ' </, is the 
CP violating phase difference between the mixing amplitude and the decay amplitude, and Ar s , the 
width difference, is predicted by the SM to be 

Ar s (SM) = 0.098 ± 0.024 ps" 1 . (4) 



The latest results from CDF with 2.8-5.2 ftr 1 luminosity are Ar s = (0.121 ±0.051) ps" 1 and $ J J^ = 
0.01 ± 0.17 [9J. Combining the results, we have Ar s = (0.134 ± 0.039) ps" 1 and $ 3 J^ = 0.21 ± 0.12, 
differing by a modest 1.7a from the SM. 

The recent like-sign dimuon charge asymmetry A^, of semileptonic 6-hadron decays from D0 is 
A b Bl = (-9.57 ± 2.51 ± 1.46) x 10~ 3 [lj for 6.1 ftr 1 luminosity. Using the SM predicted value of 
ofj = (— 4.8_ 1 ' 2 ) x 10~ 4 (we assume that the new physics contribution in the B^ system is small) and 
combining with the explicit measurement af, = —(1-7 ± 9.1) x 10 -3 [lOj . the charge asymmetry for 
"wrong-charge" semileptonic i? s -meson decay is af, = —(12.5 ± 4.8) x 10 -3 . The SM prediction is 
Og|(SM) = (2.1 ± 0.6) x 10~ 5 [3j, which is off from the measured quantity by 2.6a\ 

The three quantities Ar s , cj)f and af, are not independent. We have the following relation among 
them [TT] 



|r 12 | . , s i Ar, 

\mi z \ flm, 



a sl = M21 sin ^ = ~^~ tan ^' • ( 5 ) 



The mass difference is measured very precisely: Am s = 17.78 ± 0.12 ps . Using the central value of 
Am s and combining the measured values of a s sl and Ar s , we find a good fit with 

Ar s = 0.134 ±0.031 ps" 1 , tan <j)f = -1.66 ±0.64, (6) 

which are off from the SM predictions by 0.9<t and 2.6<7, respectively. As we are anticipating new 
physics contributions to T 12 which do not necessarily contribute to B s — > J ftp $ decays we distinguish 
the fit to fis^ from the fit to (f>f. 

In constrast, the measured mass difference in the B s system and the ratio of the mass differences 
between the B s and Bd systems can be used to extract CKM matrix elements which are in fair agree- 
ment with those extracted from other observables, although there is room for an O(20%) contribution 
to either mass difference from new physics [12^ IT3] . 

New physics can in principle contribute to both T 12 and m 12 . It is not possible to get a good fit to 
the dimuon asymmetry obtained by D0 in terms of a new contribution to ml 2 alone |14} [T5] . When 
the SM value of T 12 is used in Eq. (I5J), fitting the leptonic asymmetry requires an unphysical value for 
the phase <pf, 

sin^ 1 = -2.3 ±1.3 . (7) 

A better fit is obtainable via a new contribution to T 12 , which also can better fit the modest deviations 
of T s and Ar s from the SM. A general feature of models with new contributions to T 12 is that, in 
contrast with models which only modify m 12 , the relation cf)f = — 2/3 s does not hold |16| . 



Effective higher dimension operators offer a general approach to any short distance new physics. 
Most attempts to explain the charge asymmetry have considered new short distance contributions |17} 
HBl[I21|2ni|2Tl|221|231|221|25|. In terms of B s and B s and decays into SM light particles, one can 
incorporate new physics by writing down effective operators and finding their allowed new decay 
channels. This approach reveals the difficulty of obtaining large new contributions to T 12 , and of 
obtaining a good fit to all available data. New dimension 6 contributions to AB = 2 operators 
contribute only to ml 2 . It is theoretically straightforward to construct theories which will produce 
such contributions which are comparable to those of the SM, since the SM AB = 2 operators occur 
only at the one-loop level and furthermore are proportional to small off-diagonal Vckm elements. 
Obtaining a new short distance contribution to V 12 requires AB = 1 operators which can contribute 
to b quark decay. The SM AB = 1 operators are produced at tree level and so for new short distance 
physics to be important, there must be new tree level contributions to b quark decays which are 
comparable to the contribution of the SM weak interactions. This is difficult to reconcile with the 
many successes of the SM in predicting lifetimes and branching fractions of B hadrons. Furthermore 
any new particles contributing to b quark decays at the tree level also can produce AB = 2 operators at 
the one loop level, and so should give a large contribution to to] 2 . A recent analysis of the constraints 
on nonstandard AB = 1 operators which could contribute to I^ 2 has been given in Ref. |26j. 

It is however conceivable that new physics could also occur at longer distance scales, and would 
have escaped notice so far provided it is sufficiently weakly coupled to the SM particles. In this paper, 
we explore in a concrete model how new, weakly coupled physics at the several GeV scale can give 
a large contribution to I^ 2 while giving a contribution to ml 2 which is smaller than that of the SM. 
In particular, we assume that there is a new spin-zero particle, called £, which can mix together with 
B s and B s . Although the amount of mixing could be small (otherwise a large modification on Am s 
is anticipated), provided ( has a much larger width than the B s , AT S can be increased to match the 
experimental value even with a very small mixing between £ and B s , B s . As the amount needed 
to increase Ar s is comparable to the amount needed for the discrepancy in T s , such a model may 
explain both anomalies. Furthermore, if the mixing parameters between £ and B s , B s contain new 
CP violation phases of order unity, a large (j)f may be obtained to explain the charge asymmetry. 

This new light scalar particle £, which has a weak coupling to the SM particles, has many pos- 
sible origins from a model building point of view. In this paper, we take a purely phenomenological 
approach. We treat this particle as a generic pseudoscalar, which for instance could be a Pseudo 
Nambu-Goldstone Boson (PNGB). For example, it may behave as a familon |27| from some spon- 
taneously broken approximate family symmetry. Or it may be a meson in a hidden sector, perhaps 



composite. The £ particle may decay directly into SM particles or into other hidden states, which 
then decay into SM particles. 

To have new contributions to AT S , one should have £ mixed with B s and B s . Such mixing can 
also modify Am s . Simply from perturbation theory, one can estimate that the modifications on the 
widths are proportional to the square of the mixing angles, while the modifications on the masses have 
an additional factor proportional to the mass difference of £ and B s or B s . So, without doing detailed 
calculations, if the anomalies can be explained by mixing with another state, the contribution to Am s 
can be reduced provided the new particle mass is close to the average B s meson mass: Mq ~ m^ . 

Our paper is organized as follows. We will first describe the interactions of this new scalar field in 
Section [2j then we will diagonalize this three-particle system in Section [3j In Section [4j we perform 
a x 2 based analysis to determine the best-fit region of the model parameter space. After that, we 
discuss various viable decay channels and conclude in Section [5] 

2 Interactions of this New Scalar Field 

In this section, we will focus on flavor changing interactions of the new spinless particle with b and s 
quarks and leave its interactions with other particles for Section [5] 

For our analysis, it is convenient to take £ to interact with SM fermions dominantly through 
derivative couplings, as would be the case for a PNGB. General flavor changing interactions may be 
written as 






C = Md»C- xM c 2 C 2 + -^C^7 M (^' + si'75)^ + />.c. + 



2 fyC0"C - 2 M ( C F C ^ K (m ' " mj) + 9 ^ mi + mj) 75 J ^ + h '°' + " ' ' (8) 

Here, Mq is the PNGB mass; ipj denotes mass eigenstate SM fermions; and F is a parameter which 
could be the spontaneous symmetry breaking scale of some global symmetry. The flavor-dependent 
couplings g l y and g 1 ^ are in general complex numbers. Other couplings could also exist, but will not 
be relevant for this part of our analysis. We will consider some other interactions in Section [4] and [5] 
The general interaction terms to describe the off-diagonal couplings with second and third generation 
quarks are 

-jd^b^^ + g'ii^s-jd^t^ig^ + tii^c + h.c. (9) 

In a model where Q is related to the breaking of global flavor symmetries, we would anticipate gy ' A ~ 
gy A , if the up- type quarks and down- type quarks transform similarly. In principle, this new particle 



£ can also couple to the first-generation quarks. We assume that such couplings are small and neglect 
them. 

For the first operator in Eq. ([9]), we integrate it by parts and use the following matrix element 



dr(0\b>f<K8(p)\B a (p))=f Bt m 2 Ba , 

yielding mass mixing terms between £ and B s , B s 

e ia f 2 CB s + e- ia f(B s , 

with a = arg(^ s ) and f 2 = \g b ^ \ f Bs m 2 B jF. With f Bs « 0.231 ± 0.015 GeV [28J and m Bs 
5.3663 ± 0.0006 GeV, we have 



(10) 



fill 



/ = 0.0026 x 



6s I 



F/\g b A 

10 6 GeV 



-1/2 



GeV. 



(12) 



We work in a basis where m 12 is real in order to give physical, reparameterization invariant, meaning 
to the phase a. If a is not zero, a new source of CP-violation enters the B s and B s system. 

The C field may decay into other light particles in its own hidden sector, or into SM particles. At 
this moment, we will simply assume it has a non- negligible width T^ and come back its decays later. 
So, in the model we are considering, there are four parameters needed to compute the effects of mixing 
with £ on the B s system: Mq, T^, f and a. 

3 Diagonalization of the Mass Matrix 

The mass-squared matrix can be written in the basis (B S ,B S ,() as 



M z 



( 



\ 



m B a 

Am Bs m Bs 

e -ia p 



Am Bs m Bs e ta f 



2 \ 



in 



e ia f 2 



e -ia p 



M 2 



(13) 



C / 



which can be diagonalized by the following unitary matrix as UM 2 W = diag{m 2 ,m 2 ,,m'^}, 

/ e^l2 _ e-' 9 ^ iff \ 

V2 v^ 13 

JJ — e' fl 12 e~' fl 12 n 

V2 V2 23 

e w i2(i9 13 -e 23 ) e -^i2(-i6>i 3 -e 23 ) 1 

\ a/2 V2 L I 



(14) 



Here, the three rotation angles are 



'12 



'13 



arctan 



rn 



B s 





— sin (2a) 






cos(2a) - 

a + e 12 ) 

-Ml ' 


Am Bs m Bs (M ( 2 - 
f 4 


- m % s ) 
#23 


_ V2f 2 cos(a + e 12 ) 
m% s - M% 



(15) 
(16) 



M?\ and only kept the leading terms. From the above 



where we have assumed that / <C \m^ — m* 
equations, one can see that when / = 0, all three mixing angles are zero and the C, is decoupled from 

\ | and a = 0(1), the mixing angle 612 is of 



the B s and B s system. When / 4 ~ Ams s ujb s \M? 
order unity. While for f 4 S> Am,£ s ra _B s \M 2 



m 



-a and #13 ~ 0. 



C -™b.\> V12 
In the diagonalized basis {B\, B2, B 3 ), the three mass eigenvalues are 



mi 



TO2 



m 3 



m Bs 



m Bs + 



Amg s 
2 

Am Bt 



COS (2^12) 
COS (2^12) 



& (M| 



m 



B. s 



? 2 
7 23 



2m Bs 

{Ml 



m 



Bs 



2m jBs 



/ 4 
Mc+ M c (M|-m| s ; 



(17) 



Including the decay width matrix and working in the basis (B\,B2,B 3 ), we need to diagonalize the 
following effective Hamiltonian 



H 



( 



\ 



Vfl\ 



nri2 



\ 



m 3 J 



U 



( 



\ 



-r 



12 







-V s 

1 12 

r, 









\ 







t/t. 



Considering the relative phase between M± 2 and Tf 2 is small in the SM, we neglect the phase of 
Tf 2 from now on. We also choose the phase convention such that Amg s and Tf 2 are both positive 
quantities. Noting that r s ,rf 2 <C AniB 3 , we use the ordinary perturbation theory to calculate the 
eigenvalues and eigenvectors. The three eigenvalues are calculated to be 

\2' 



VB S , L 
VB StH 



mi + 



(i 



127 



1712 



m 3 



O 

i 

2 



\Am Bs 



s \2 
VI) 



(t 



\Am,B s 



[F s + T c 9 2 13 + Ti2 cos(20i 2 )] 

[r s + r c 0| 3 - r 12 cos(2^ 12 )] 



[r< 



r c (^ 3 + ^ 



2?J 



(19) 



Neglecting terms suppressed by \/{m\ — m.3) and l/(m 2 — m 3 ), the formulae of the mass eigenstates 
of the two lighter one in terms of flavor eigenstates are 

B _ e i0 ^ ( fl 13 g 23 r c -sin(2fl 12 )r 12 \ e- i0 ™ ( giagggTg - sin(2g 12 )r 12 \ 

S ' L V2 V 2(mi-m 2 ) / s ^2 V 2(mi - m 2 ) / s 

+ 0(013,023) C, 

R e ie ™ ^ , 0i3g2 3 r c -sin(2g 12 )r 12 ^ e~^ / _ %0 23 r c - sin(2g 12 )r 12 ^ 

s '" " V2 V 2(mi -ma) J * + V2 \ 2(mi - m 2 ) J ^ s 

+ O(0i3,0 23 )C, (20) 



The mass and width differences of the heavy state and the light state are 

_ (0|s - #13) ( M c ~ m IJ 

Am s = mg -m L = cos (26» i2 ) Am Bs , (21) 

2 m Bs 

Ar s = r L - v H = r c (0? 3 - e 2 23 ) + 2 r i2 cos(2#i 2 ) . (22) 

The average width is 



A 2 1 (ft 

r s = r s + ^±^r c . (23) 



From Eqs. (21) and (22), we can see that compared to the width difference the mass difference can 
be suppressed by an extra factor [Mc, — mB s )/mB a if M^ is close to rns s . This fact provides us the 
possibility of increasing the width difference without changing the mass difference too much. 

Neglecting the effects of £ in the B s and B s oscillation, we use the traditional formula to calculate 
the charge asymmetry 



T(B S -)- n+X) - T(B S -)■ fi-X) 


|E|2 _ 1 <?|2 
Igl Ipl 


T(B S -> n+X) + r(5 s -)■ /i-X) 
2 0i 3 23 r c 2sin(20i 2 )ri 2 


1212 _i_ 1912 
1 gl "I" Ipl 


mi — m 2 mi — m 2 
2%0 23 r c 2sin(2# 12 )r 12 





(24) 
Am s Am s 

When a = and no new CP violating physics exists, we have Qyi = #13 = and hence af, = 0. In 
order to compare to the observables both in the charge asymmetry and in the B s — > J/ip $ decay, we 
calculate the phase cj)f and f3 s in our model 

sl _ 2g 13 g 23 r c -2sin(2g 12 )r 12 

0S " r c (^-^3) + 2ri2co8(2e 12 )' [b) 



4 Fit to Observables 

In this section, we want to use this new model to fit the five observables in the neutral B s meson 
system. We first summarize the various experimental values and the SM predictions in Table [T] 





Experimental 


SM prediction 


Am s 


(17.78 ±0.12) ps- 1 


(19.6 ±2.2) ps" 1 


AT S 


0.134 ±0.031 ps" 1 


(0.098 ± 0.024) ps" 1 


r s 


0.680 ±0.012 ps" 1 


(0.654 ±0.008) ps" 1 


tan 4>f 


-1.66 ±0.64 


0.0042 ± 0.0014 


$** 


0.21 ±0.12 


0.018 ±0.001 



Table 1: The experimental values and SM predictions for the five observables considered in this paper. 
We have four model parameters, Ms, f, a, IV, to fit the five observables, Am s , Ar s , F s , tan^ 1 



and j3 s .To quantify the goodness of fit from the new physics, we define the following x 2 



X' 



E 



(Of 



model 



QCXP\2 



a, 



SM 



± (T. 



(27) 



cxp 



with Oi represents the five observables. Neglecting new physics contributions or setting / = 0, we 
have x 2 (SM) = 14.0, which indicates a large discrepancy between the SM and those five observables. 

We first look at the case with an approximately massless PNGB. Since the two-body decay process 
b — > s + C from the spectator model is open, we should also consider constraints from the decay width 
of Bd- The experimental measured value is Ts d = 0.634 ± 0.004 ps . However, the SM prediction 
for this quantity has a large uncertainty. For example, the decay constant fs d = 190 ± 13 MeV [28J 
from the lattice QCD calculation, which gives around 14% uncertainty to r# d . To be conservative, 
we neglect other possible uncertainties and require the new physics contribution to r# d to be less 
than 0.09 ps . Other B meson and B hadrons do not constrain the parameter space further, because 
the relative experimental errors of their widths are higher than T# d and the theoretical errors are 
comparable. 

For M^ < rriB d — mriK with tjik as the kaon mass, the two-body decay width of this channel is 



calculated to be 



ml\q b n 2 (m 2 B -M 2 ) 2 m lf (ml - M 2 ) 2 



16ttF 2 m? B 16 7r/o mi 






771, 



Neglecting M^ and requiring r spec < <5Ts d = 0.09 ps _1 , we have a constraint on / < 1.1 x 10~ 3 GeV. 
Although a pretty good fit can be found for the approximately massless case, the allowed region is 
ruled out by the constraint from rs d . The best fit has x 2 = 2.0, and we present 68% and 90% contours 
around the best-fit region of our model parameters in Fig. [TJ We therefore conclude that the mass of 
the C field cannot be light compared with the 6-quark mass. 
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90% 



O 



o.ooo L -1 00 ° r 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 1: Left panel: the best-fit region in the / and a space for a light £ mass and a fixed width. The 
global minimum has \ 2 = 2-0. The two contours have 68% and 90% C.L., respectively. The region 
above the black solid line is ruled out due to decay width of B^. Right panel: the same as the left 
panel but in the T^ and a plane. 



For Mq > ms d — rriK, the two-body decay channels are forbidden, but there are still three-body 
decay channels open. The three-body decay width is related to the width T^ because an off-shell £ 
mediates the three-body decay. To be concrete, we assume that £ can decay into two light scalar 
fields a via the interaction k£o 2 /2, where k has mass dimension one. The light scalar field a can 
subsequently decay back into SM particles. Detailed discussions about possible decay channels will be 



presented in Section [5] We have the two-body decay width T^ as 



T c (C^2a) 



(29) 



32 it M c 

Defining the Yukawa coupling A = m& \g h ^\jF and neglecting the mass of a, the three-body decay 
width, Y^iBd — >■ K + a + a) through an off-shell £, is calculated as 

A 2 M 2 



r. 



32 ?H 



S(M^,m Bd ,mK) ■ 



(30) 



For A = 10 5 and \k\ = 1 GeV, we calculate this decay width using Calchep [29] and show it in Fig. 
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Figure 2: The three-body decay width of Bd — >• K + o + a for A = 10 and |/s| = 1 GeV. 
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Using Eq. (29) and the relations / = \g^\ Ib s itIb /-F 1 and A = rrn, \g°?\/F, we have the three-body 



decay width as a function of / and IV 



/ 4 m 2 M c r c 



S(M c ,mB d ,m K ) 



(31) 



In the following, we will impose the constraint T3 < <5r# d = 0.09 ps _1 . 

We present the allowed parameter space in Fig. [3} In the left panel, we first fix the decay width 
as r,j = 0.001 GeV, and then calculate the total x 2 by marginalizing a. We have found a region 
of parameter space providing a much better fit to the five quantities in Table [I] than from the SM. 
The best fit can have x 2 = 2.0. The right panel is similar to the left one but with T^ = 0.01 GeV. 
As can be seen from those plots, the £ is preferred to have a mass close to the B s meson. This is 
because when M^ ~ rriB s , one can have fairly large changes to the quantities Ar s and aj without a 



10 



large contribution to Am s . Comparing those two plots, one can see that the plot with a larger T^ has 
more parameter space ruled out by the three-body decay width, which can be understood simply from 



Eq. (|31|). To illustrate the goodness of our fit, we report the results for one point of our parameter 

(32) 



space 

M c = 5.2GeV, / = 0.0023 GeV, T c = 0.0025 GeV , a = 1.10. 

For those numbers, we have the following model prediction 



Am™ d = 17.23 ps" 1 



Ar s nod = 0.125 ps- 1 



=rmod 



0.690 ps~ 



tan 



i si mod 



-0.70, 



pj/1* 



mod 



0.13, 



(33) 



which has good agreement with the experimental measured values and has a total x 2 = 3.2 compared 

-5.0 x 10~ 3 for this point. 



to x = 14.0 in the SM. We also report the charge symmetry as aj 
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Figure 3: Left panel: the best-fit region in the M^ and / space for a fixed width T^ = 0.001 GeV. The 
orange contour has 68% C.L. after minimizing \ 2 hi terms of a. The best fit has x 2 = 2-0. The gray 
region is ruled out by the two-body decay width of B^ when M^ < niB d — mx- Three-body decays 
do not rule out the best-fit region. Right panel: the same as the left panel but for T^ = 0.01 GeV. 
The best fit has \ 2 = 5.4. The blue region is excluded by requiring the three-body decay width to be 
below the error of 6Tg d = 0.09 ps . 



We present the best-fit region in the T^ and / plane by fixing a specific £ mass M^ = 5.2 GeV in 
Fig. |4| In the left panel of this figure, we still treat a as a floating parameter. The best-fit region 
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Figure 4: Left panel: the best-fit region in the T^ and / space for a fixed mass M^ = 5.2 GeV. The 
orange contour has 68% C.L. after minimizing \ 2 m terms of a. The best fit has \ 2 = 2.0. The blue 
region is excluded by requiring the three-body decay width below the error of 6Ts d = 0.09 ps . Right 
panel: the best-fit region in the a and / space for Mq = 5.2 GeV and T^ = 0.002 GeV. The best fit 
has x 2 = 2.3. 



prefers T^ within 10 -4 — 10 -2 GeV. The best-fit region is not ruled out by the three-body decay width 
of Bd- In this plot and to have weak coupling between Q and its decay products, we don't extend the 
width of C to be above 1 GeV, which is around 20% of its mass. In the right panel of this figure, we 
present the allowed region in a and / for fixed values of Mq = 5.2 GeV and T^ = 0.002 GeV. 

5 Discussion of New Meson Decay Modes and Conclusions 

From the best-fit region in the left panel of Fig. |4j the width of this new pseudoscalar particle should 
be above 10 -4 GeV. We have only considered the total width constraint on various B meson decays 
so far. For some decay products from £, more stringent bounds may be applied. Again, considering 
only two-body decays of Q to SM particles and neglecting the final state masses, the coupling A( of Q 
to SM particles should be above ~ 0.1. 

We first consider the situation where the C, directly decays into two SM particles. Due to the strong 
constraints on branching ratios of B^ to K° plus e + e~, / u + ^ , uu, 7r + 7r~, K~tt + and K + K~ , we are left 
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with options £ — > e + r , /i + t , t + t , cu, and cc. The anomalous magnetic dipole moments (g£ — 2)/2, 
which is approximately \ 2 -/(4'7T 2 )(m 2 /M?) 2 , are constrained to be smaller than ~ 10 , ~ 10~ 9 , and 



fsJ 



10~ 2 for the electron, muon and tau, respectively. So, the two channels (, — > e + r~, fi + T" are ruled 
out. For the t + t~ channel, which is also pointed out recently in Ref. |30j, one may worry about 
direct searches, for example at LEP II. While the production cross section e + e~ —> t + t~Q — > 2t + 2t~ 
with M^ = 5.0 GeV and A^ = 0.1 is ~ 1 fb and the total luminosity at LEP II is around 0.7 fb . 
One may also worry about the modification on the width of the Z boson decays. The decay width of 
T(Z — > t + t~Q) « 0.013 MeV for Mq = 5.0 GeV and A^ = 0.1, which is also below the measurement 
error of mz = 91.1876 ± 0.0021 GeV. Therefore, we conclude that the t + t~ channel is allowed. 

The AC = 1 couplings Qcu can contribute to an effective AC = 2 operator which is highly 
constrained by D — D° mixing, giving rise to a mass difference of order 



Am D ~ X 2 c (f 2 D m D )/Ml (34) 

compared with the measured value of |31j 

Am D = (1.6 ± 0.5) x 10~ 14 GeV , (35) 

which constrains A^ to be less than ~ 2 x 10~ 6 . We also note that if the coupling of C c u is CP- violating, 
the constraint from D° and D mixing is even more constraining. 

For the CC75C couplings, we have mixing of £ with the rj c which gives rise to Sm rtc ~ 0.3 MeV, 
which is below the experimental error m Vc = 2980.5 ± 1.2 MeV. After fragmentation, we anticipate 
C -»• DD{n's). From spectator decays, Br(B s ->■ D~D+), Br(B d -+ D~D+) and Br(B s -»■ J/tpr}) 
are estimated to be around 10 -4 in our model, which is below but close to the current experimental 
errors [32J, (33J [HI] . So, we have two allowed two-body decay channels of (: t + t~ and cc. 

Also possible are three body decays, e.g. £ ~~ * cua, with a a new stable particle. If a differs from 
its own antiparticle then it is possible to avoid an excessive contribution to AC = 2 operators. We 
anticipate £ - > D(ir's)X with X denoting missing particles. Making a three-body decay consistent 
with a large enough T^ is difficult however, due to phase space constraints. 

A third possibility is to have £ first decay into other lighter particles in the hidden sector and then 
have those lighter particles decay back to the SM particles. If £ decays into two identical particles a, 
more possible channels are open for a decay. The a field can decay into e + r~ and /j> + t~ if its couplings 
A a ^ 10~ 4 to evade the constraints from the anomalous magnetic dipole moment. The ct of the o 
field is then estimated to be above ~ 1 \im. Although the modes to t + t~ and cc are kinematically 
forbidden because we need M a < M^/2, the mode cu is allowed. However for such a decay to take 
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place without a significant displaced vertex requires a coupling which is too large to be consistent 
with D° — D mixing, unless a is not its own antiparticle. The decays a — > 7r + 7r~ or K + K~ are 
also allowed because the exclusive searches for decay products of Bj have not yet covered this kind 
of high-multiplicity final states like Bj — > K°7r + 'K~7r + 'K~ . Decays into K~ it + are allowed, but if a is 
its own antiparticle then avoiding a AS = 2 contribution to K° — K mixing will imply a significant 
displaced vertex. 

We also have the option of (, decay into two different light particles a\ and 02- Now, we can have 
more combinations of final state particles with a\ and 02 decaying into different SM model particles. 
More specially, we can have one particle a<i to be semi-stable and missing particle if its couplings to 
SM particles are weak enough to escape the detector. However, we can not allow both a\ and ai to 
be (semi)stable particles because of the constraint from Br(Bd — > K®vv) < 1.6 x 10 . 

A fifth possibility, which may be difficult to constrain, is that £ decays into 2a\ followed by decays 
of the a\ particle into SM particles and a hidden particle 02, where 02 may be semistable and escape 
the detector. 

We summarize the simplest allowed decay modes of £ in Table [2] The new decay modes should 



Direct decay 



(^2a 



C — > a\ + a 2 



Decay Modes 
t+t-, D D{rfs), D(tt's)X 



2T + 2e~, 2t+2/i-, 2D+271-, 2tt + 2tt-, 2tt + 2tt-, 2K~2tt + , 2K+2K~ 



X + (r+e~, t + /jT, D + Tr~, 7r + 7r~, 7r+7r~, K^tt + , K + K~ 



Table 2: Some allowed decay modes of £, with X representing missing particles. Other possible final 
states with different combinations of charges are also allowed. For example in the Q — > a\ + a% case, 
different combinations of final states in the parenthesis are also allowed and not shown here. 

account for approximately 3.5% of the total width of the B s using the best-fit region in Fig [3| So, if 
future experimental results find those new decay modes of B s in Table[2]but not for Bj, our prediction 
for this new state £ would be confirmed. 

Before we end our paper, we make comments about the second operator in Eq. ([9]). This operator 
can induce flavor changing neutral current type decays of the top quark: t —> c + C,. Depending on the 
final states of C, decay, there may exist a decay channel like t — > c + t + t~ , with the invariant mass of 
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the t pair around 5 GeV. The branching ratio of this new decay channel is calculated as 

From Eq. (12) and Fig. bl we see that if g^ ~ g h j[ the explanation of anomalies in the B s symmetry 
in the familon framework predicts this branching ratio should be around 10~ 7 , which may be tested 
at the LHC. 

In summary, we have shown that several anomalies in the B s system can be explained simultane- 
ously if there is a new pseudoscalar, the £, mixing with B s mesons. From fitting to the five observables: 
Am s , Ar s , T s , (f)f and j3 s , this pseudoscalar is predicted to have a mass around 5 GeV and a width 
around 10 -3 GeV. Many viable decay modes can be found for the new pseudoscalar. As a result, the 
Bfa and B^ should have non-negligible branching ratios of around 10 -4 into a kaon (or other fragmen- 
tation modes of sd and su) plus the decay modes of the (, such as those listed in Table Likewise 
the B s should have a branching fraction of around 10 -4 into ss fragmentation states such as <fi or r], 
plus the decay modes of the £. In addition, the B s branching fraction into the £ decay modes (with no 
additional particles) should be about 0.035. Our model can be tested at the Tevatron and the LHC, 
particularly LHCb. We also emphasize that <pf ^ —2(3 S in our model, in constrast with models 
in which new physics contributed only to ml 2 . More precise measurements of those two quantities 
would therefore distinguish our explanation of anomalies in the B s system from other approaches. Our 
motivation for this model is to better fit experimental data, and we have not attempted to justify the 
naturalness of the model or to discuss the theoretical implications of such a new particle or particles, 
other than to make sure the model can be consistent with other experimental results. However we note 
that light pseudoscalars resulting from spontaneously broken approximate symmetries are common in 
models of flavor physics and/or dynamical electroweak symmetry breaking at the weak scale, as well 
as in models with new strongly coupled hidden sectors, which have a variety of theoretical motivations. 
It may seem somewhat surprising or coincidental that an exotic hidden meson should have a mass 
which is rather close to that of the B s , however this is the simplest viable way we are aware of to get 
a large contribution to V 12 . 
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